A Twisted Ladder: relating the Fe superconductors to the high Tc cuprates 



o 
o 



a ■ 
o ■ 

u • 

^: 
^ : 

a ■ 



I 

O 

o 



> 

o\ 
o 

o 

On 
O 



X 





(a) 



kx/a 
(b) 



E. Berg^, S. A. Kivelson^, and D. J. Scalapino^ 
^Department of Physics, Stanford University, Stanford, CA 94305-4045, USA 
^Department of Physics, University of Cahfornia, Santa Barbara, CA 93106-9530, USA 

(Dated: May 7, 2009) 

We construct a 2-leg ladder model of an Fe-pnictide superconductor and discuss its properties and 
relationship with the familiar 2-leg cuprate model. Our results suggest that the underlying pairing 
mechanism for the Fe-pnictide superconductors is similar to that for the cuprates. 



An important question has been raised by the discov- 
ery of high temperature superconductivity (HTS) in the 
Fe-pnictides: Is there a single general mechanism of HTS 
which operates, albeit with material specific differences, 
in both the Fe-pnictides and the cuprates (and possi- 
bly other novel superconductors) , or do the cuprates and 
the Fe-pnictides embody two of possibly many differ- 
ent mechanisms of HTS? This question is complicated 
by one of the perennial issue of the field: To what ex- 
tent is it possible to understand the properties of the 
strongly correlated electron fluid in the Fe-pnictides, the 
cuprates, and other materials from a weak or strong 
coupling perspective, given that the materials exhibit 
some features that appear more natural in one limit and 
some that are suggestive of the other. The Fe-pnictides, 
like the cuprates, are likely in the intermediate coupling 
regime [E S B 3 j which is difficult to treat theoretically. 

In the present paper, we introduce a model of a two-leg 
ladder with an electronic structure chosen to reproduce 
particular momentum cuts through the Fe-pnictides band 
structure. We study the magnetic and superconduct- 
ing properties of this model numerically using the den- 
sity matrix renormalization group (DMRG) '5j method 
which allows us to treat the "intermediate coupling" 
problem essentially exactly. Although the existence of 
short range, spin gapped, antiferromagnetic correlations 
in the half-filled two-leg Hubbard ladder and the power- 
law d-wave-like pairing correlations in the doped ladder 
were initially unexpected this behavior is now well 
understood 0, H 0| ■ Thus it is not surprising that a two- 
leg ladder model of the Fe-superconductors exhibits short 
range antiferromagnetic correlations and power-law pair- 
ing. However, it is interesting to see how such a model 
captures particular aspects of the physical properties of 
the Fe superconducting materials [I0|, specifically their 
spin and pairing correlations. 

We find that, in the limit of zero temperature, this 
model has a diverging superconducting susceptibility 
with a pair structure of a form which is the quasi-one 
dimensional version of the structures which have been 
proposed on the basis of both weak 11 , I^, 13. 14. 15. ijj 
and strong fl7'| coupling calculations for the 2D system. 
We also find that the dominant magnetic correlations are 
"stripe-like" with wave vector (0,7r), reminiscent of the 
magnetic structure that is seen in the undoped parent 



FIG. 1: (a) Heisenberg model showing local Fe spins coupled 
by comparable nearest neighbor Ji and next nearest neighbor 
J2 exchange couplings, (b) 5-orbital tight binding Fermi sur- 
faces with the main orbital contributions shown by the colors: 
dxz (red), dyz (green), and d^cy (yellow). The arrows illustrate 
the type of scattering processes which give rise to pairing in 
the fluctuation exchange calculations. 



compounds (and sometimes coexisting with superconduc- 
tivity) in the Fe-pnictide superconductors. In addition, 
we show that the Fe-pnictide ladder is a "twisted" version 
of the two-leg Hubbard ladder. Upon untwisting the lad- 
der, a correspondence is found with the familiar results 
from cuprate related studies: the superconductivity in 
the Fe ladder corresponds to "d-wave-like" correlations 
in the cuprate Hubbard ladder, and the stripe-like mag- 
netic correlations are transformed into antiferromagnetic 
correlations with wave vector (tt, tt). This supports the 
idea that there is a single, unified mechanism at work in 
these two families of HTSs. 

As in the case of the cuprates, both strong coupling and 
weak coupling models have been proposed to account for 
the magnetic, structural and superconducting properties 
of the Fe-pnictides. From a strong couphng perspective, 
the observed magnetic and structural phase transitions 
and much of the dynamical magnetic structure observed 
even in the superconducting phase are thought of as aris- 
ing from frustrated quantum magnetism [la. [l9j . Here, 
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as shown in Fig. [T^, one has for the undoped system 
a Heisenberg model in which the second-neighbor anti- 
ferromagnetic interaction, J2, is comparable to or larger 
than the nearest-neighbor exchange interaction, Ji. For 
the doped system, one has a. t — Ji — J2 model. Then, 
as discussed in Ref. [13], a mean- field analysis shows 
that the dominant pairing channel involves intra-orbital 
{dxz,'\, dxz,i) and {dyz,'\, dyz,\) pairs, combined to form an 
Aig superposition. A similar pairing structure was found 
in Ref. [20] in the large U limit of an exact diagonaliza- 
tion study of a a/8 x -\/8 cluster. While both of these 
studies only took into account the dxz and dyz orbitals, 
their finding concerning the importance of intra-orbital 
pairing is relevant to the model we will discuss. 

Alternatively, the weak coupling picture begins by con- 
sidering the energy bands and the Fermi surface. Fig. 
[T}) shows the Fermi surfaces for a five-orbital tight bind- 
ing fit I id] of a density function theory (DFT) calcula- 
tion [2l| , neglecting the effects of dispersion in the third 
direction, perpendicular to the Fe-pnictidc planes. Here 
the color indicates the relative weights of the 3d {dxz^ 
dyz, dxy) orbitals which form the main contribution to 
the Bloch wavefunctions on the Fermi surfaces. The dxz 
and dyz orbitals provide the dominant weight on the hole 
Fermi surfaces ai and 02 around the F point, while on 
the j3i and (32 electron Fermi surfaces one has (^dyz, dxy) 
and {dxz, dxy) contributions, respectively. 

From a weak coupling perspective, it is the nesting 
properties of the a and P Fermi surfaces that give rise 
to the spin density wave (SDW) instability which ac- 
counts for the magnetism of the undoped parent com- 
pounds. For the doped system, fiuctuation exchange 12l . 

i3, H m and renormalization group calcu- 
lations find that magnetic fluctuations near (tt, 0) and 
(0, tt) lead to a pairing instability. In this case, similar 
to the strong-coupling model, the pairing arises from the 
type of scattering processes illustrated in Fig. [lb. Here, 
a (k t,~k I) pair on the dxz region of the ai Fermi 
surface is scattered by a Q = (0, tt) spin fiuctuation to 
a (Q + k' I, Q — k' |) pair on the dxz part of the (32 
Fermi surface. Similar (tt, 0) processes involving dyz pair 
scattering occur between ai and /3i. 

Thus, both strong-coupling and weak-coupling ap- 
proaches have been used to describe these materials. 
However, what is needed is an approach which allows 
one to treat the intermediate coupling regime. Here, us- 
ing DMRG, we address this issue for a caricature of the 
original problem which focuses on the dxz orbital ai — (32 
scattering process for ky = and ky — n states near 
the Fermi surface. These scattering processes can be de- 
scribed by the two-leg ladder shown in the inset of Fig. [2] 
with a Hamiltonian given by 
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FIG. 2: Band structure e{kx,ky) of the ladder with t\ = 
—0.32, t2 = 1 and ts = —0.57. The solid and dashed curves 
correspond to ky = and tt, respectively. The black dotted 
line corresponds to the Fermi energy for a filling of one elec- 
tron per dxz orbital. The dots represent the bands from a 
DFT calculation [2^ for the 2D lattice, cut through ky = 
and TT. The inset shows the ladder with the one electron hop- 
ping matrix elements ti,t2 and t^. 



-2i3 ^ dj^i^(ij+l,2cr 



h.c. + uY^ (1) 



Here, a = 1, 2 is the leg index, a =t, | is the spin index, 
there are leg ti , rung t2 and diagonal one-electron hop- 
ping matrix elements and an on-site Coulomb interaction 
U. The factors of 2 in front of t2 and ^3 take into account 
the periodic boundary conditions transverse to the lad- 
der. The hopping strengths were fitted to the DFT band 
dispersion of the ai and (^2 pockets calculated in Ref. 
[2l|, cut through ky = and ky = tt near the Fermi sur- 
faces. At these points in momentum space, the Bloch 
wavefunctions have a dxz character. The parameters we 
used were ti = —0.32 and = —0.57, measured in en- 
ergy units in which t2 — 1- The ladder bandstructure is 
shown in Fig. [21 where it is compared to the DFT disper- 
sions near the Fermi energy. The dotted line marks the 
location of the chemical potential for the half-filled, one 
electron per dxz orbital case. To capture the intermedi- 
ate coupling aspect of the physics, we will take the onsite 
Coulomb interaction [7 = 3. This value of the interaction 
is smaller than the overall bandwidth of the ladder, but 
larger than the Fermi energies of the hole and electron 
pockets relative to their values at /cj, = 0. 

DMRG calculations were carried out for a 32 x 2 sys- 
tem. As expected, there is a spin gap, and the doped 
system exhibits power-law pairing correlations. At half- 
filling, in the presence of an externally applied magnetic 
field on the first site of the lower leg, we find the striped 
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FIG. 3: (a) Spin structure for the undoped ladder, for U /t2 = 
3. The lengths of the arrows indicate the measured values of 
{S^{lx)) in a calculation in which an external magnetic field 
has been applied to the lower leftmost site. (S'^ilx)) decays 
exponentially with a correlation length of about four sites, 
(b) Fairfield structure for the (n) = 0.94 doped ladder. Here, 
a pairfield boundary condition has been applied to a rung on 
one end of the ladder and the induced rung pairfield {Pi^) 
versus Ix is shown on a log- log plot. The amplitude of the 
singlet pairfield across a rung and along a diagonal at Ix ~ 10 
are shown in the inset. The error bars in the pair fields were 
estimated from the extrapolation of the DMRG truncation 
error to zero. The pairfield amplitude along a leg is less than 
10"^ 



spin pattern shown in Fig. [3]a,. The spin correlations de- 
cay exponentially with a correlation length of about four 
sites and a spin gap A^. = 0.14. 

For the lightly hole doped case, (n) = 0.94 (4 holes), 
we find a spin gap Ag = 0.07 and power-law pairfield 
correlations. Here, a pairfield boundary term 



with Ai — 0.5 and 



P 



'l,2i 



h.C. 



*l,li"l,2T 



(2) 



(3) 



was added to the Hamiltonian. This acts as a proximity 
coupling to a rung on the end of the ladder. 

The two leg ladder is in a phase with a single gapless 
charge mode, and the long range correlations are power 
laws characterized by the Luttinger parameter Kc- For 
Kc > J the superconducting susceptibility is divergent. 



and for Kc > | it is dominant over the charge density 
wave susceptibility. In addition, if Kc > \ , the boundary 
pairing term of Eq. ^ is relevant under a renormaliza- 
tion group flow [2j, |25| and the effective boundary con- 
dition is described by perfect Andreev reflection. In this 

1 

case, the induced pair field decays a s [ L tan (fj^) | "^'^ , 
where L is the length of the system [20| . For the stated 
parameters, the induced expectation value of the rung 
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pair amplitude (P^fe)) - \"';,,iT"'ix,2j, - "'i,,ii"'!,,2T 
was measured throughout the system, and is shown in 
Fig. on a logarithmic scale. The slope of the curve 
\og{{P\lx))) VS. \og{lx) gives Kc — 0.5, so the ladder is 
near the border of the superconducting [Kc > 0.5) phase. 
This conclusion is supported by calculations with longer 
systems (up to 64 x 2). For electron doping, we found 
similar results with slightly smaller values of Kc- 

The inset of Fig. shows the amplitude for removing 
a singlet pair from two sites at a position 10 sites away 
from the boundary. In order to interpret these results, 
we perform a BCS mean field treatment of the ladder. In 
such a treatment, the amplitude for removing a singlet 
pair from two sites separated by {Ix, ly) is 



A{lx,ly) = {di+i^^j+i^^di^ji- di+i^j+i^idi^j^) 



1 vA(k) ,^.1 

£;(k) ■ 



(4) 



Here £; (k) = [e (k) - /z] + | A (k) P, iV is the number 
of sites, and the prime on the k sum implies that it is 
cut off when |e(k) — /i| > wqi where ojq is a cutoff energy. 
We assume that the gap function has a simple form in 
which A(k) is approximately constant for k in the neigh- 
borhood of a given Fermi point, but with a sign change 
between the ky — Q and tt bands. A reasonable estimate 
of the relative magnitudes of the gaps on the two bands 
is obtained by requiring that, due to the strong onsite 
Coulomb interaction, the onsite pair amplitude vanishes: 



A (0, 0) ~ 2 ^ N (0, ky) A [ky) In 



A{ky) 



0. (5) 



Here N{0,ky) and A{ky) are the density of states at 
the Fermi energy and the gap function, respectively, 
of the ky ~ and tt bands. We have assumed that 
A{ky) ^ luq ^ W where W is the bandwidth. From 
Eqs. (jll and ([5|), it follows then that the ratio of 
the amplitude for removing a singlet pair from a diag- 
onal bond to removing it from a rung, A (1, 1) /A (0, 1), 
is approximately i (cos fc^ (0) + cos A:_f (tt)) « 0.9. Here 
kp (0) and kp (tt) are the kx Fermi momenta for ky — 
and TT, respectively. Similarly, the ratio of the 
leg to rung amplitude A (1, 0) /A (0, 1) is approximately 
^ (cos kp (0) — cos kp (tt)) « —0.02. These amplitude ra- 
tios are in overall agreement with the numerical results 
shown in the inset of Fig. [A (1, 0) is less than 10""^, 
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FIG. 4: (a-b) Mapping of the Fe Hubbard ladder, Eq. (1) to 
the Hubbard ladder used to model the cuprates. (a) Twist 
every other rung as indicated, (b) Change the sign of the d^z 
orbital on the shaded sites, (c) The (0, n) short range spin 
correlations of the FeAs ladder (shown in Fig. [3^) become the 
familiar (vr, tt) spin correlations of the usual Hubbard ladder 
after the twists and orbital sign changes shown in (a) and (b). 
(d) The pairing correlations shown in the inset of Fig. 3b are 
mapped to the familiar cuprate ladder d^2_j,2-Iike pairfield. 



which is essentially zero within our numerical accuracy.] 
The magnitude of the pairing amplitude depends on the 
strength of the boundary proximity pair field Ai. Note 
that the rung and diagonal pairfields are in phase with 
each other. A similar pairfield pattern was found in weak 
coupling spin-fluctuation calculations (l6| . In this case, 
there were also dyz-dyz orbital pairfield correlations, ro- 



tated by 90°. The relative phase of the dx 



and 



dyz — dyz pairfields determines whether the order param- 
eter has Aig or Big symmetry. The present study can 
not address this issue since it lacks the dyz orbitals [27 1. 



At first sight the striped (0, tt) spin structure of the 
undoped ladder and the pairing amplitudes of the doped 
ladder appear different from the (7r,7r) spin and dj;2_y2- 
like correlations familiar for the ladders used to model 
the cuprates. However, a closer look at the present model 
reveals that it is in fact identical to that used to model 
the cuprates. 

The mapping between the Fe model we have been dis- 
cussing and the "usual" two-leg Hubbard ladder is il- 
lustrated in Fig. m First, every other rung is twisted 
(Fig. 2^), interchanging the two sites at the ends of the 
rung and, leading to a ladder having a hopping ^3 along 
the legs and the weaker hopping ti along the diagonals. 
Then, as illustrated in Fig. Hb, the phases of the orbitals 
on the shaded sites are changed by —1. This leads to the 
usual Hubbard model with iiog = —2^3 — 1.14, irung — 
2t2 = 2 and a next-nearest neighbor t' = —ti = 0.32. 
After the twist and orbital phase changes, one finds the 
familiar (tt, tt) spin correlations shown in Fig. H]: and the 
d2:2_j^2-like structure for the pairing correlations shown 
in Fig. Hli. 

Finally, we note that the parameters for the FeAs lad- 
der place it near a region of enhanced pairing for the 



corresponding Hubbard ladder case (8[. There, it was 
found that the pairing correlations are enhanced when 
the parameters are such that the Fermi level is located 
near the top of the bonding band and the bottom of the 
anti-bonding band (^rung/^ieg ^ 2). The proximity to 
this point also makes the pairing correlations sensitive to 
the bandstructure parameters. We have found that while 
using different bandstructure parameters and a different 
U does not change the behavior found here qualitatively 
(e.g., the existence of a spin gap and the pair structure 
is robust), the strength of the long-range pairing corre- 
lations is very sensitive. For instance, for the bandstruc- 
ture parameters used in the present paper, the pairing 
correlations are reduced upon increasing the interaction 
to U = i. For this case, we get Kc ~ 0.3. In addi- 
tion, if we use the band parameters ti — —0.77, t2 = 1, 
t3 — —0.65, obtained from a fit to a 2-orbital model [l^ . 
and U — 4.6, we get a significant enhancement of the 
pairing correlations, and Kc ~ 0.75. The sensitivity of 
the pairing correlations to the band parameters and to U 
may be an artifact of the one dimensional nature of the 
ladder [2^. We believe, however, that the basic super- 
conducting mechanism in the ladder is operative in the 
FeAs superconductors. 

In summary, this analysis began by constructing a lad- 
der model for the Fe-pnictides with hopping parameters 
chosen to fit DFT calculations. Following the results of 
the RPA spin-fluctuation calculations and the mean-field 
t — Ji — J2 and exact diagonalization studies which in- 
dicate that the dominant pairing involves (0, tt) dxz to 
dxz and (tt, 0) dyz to dyz pair scattering processes, only 
the dxz orbital was kept. Using the DMRG method, we 
found short range magnetic correlations which mimicked 
the "striped" (0, tt) magnetic order in the undoped par- 
ent system, and pairing correlations with a structure that 
agreed with both the weak coupling RPA calculations 
and the strong coupling results. The Fe-ladder turns 
out to be simply a twisted version of the usual cuprate 
Hubbard ladder, with parameters near the regime of en- 
hanced pairing. We believe that this is not an accident 
but arises from the close connection between the basic 
physics responsible for high Tc superconductivity in the 
cuprates and the Fe-pnictides. 
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